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LONG-TERM  GOALS 

The  prediction  of  water-leaving  radiance  in  coastal  waters  is  strongly  dependent  on  a  quantitative 
prediction  of  the  depth-dependent  distribution  of  Colored  Dissolved  Organic  Matter  [CDOM]  in  the 
water  column.  The  goal  of  this  project  is  to  support  and  synthesize  laboratory  and  field  experiments 
conducted  by  the  ONR  Environmental  Optics  Program.  Our  part  of  this  larger  project  is  to  develop  a 
quantitative  understanding  of  the  physical,  chemical,  optical,  and  biological  interactions  impacting 
CDOM  cycling,  and  to  codify  this  understanding  into  a  numerical  simulation,  EcoSim  2.0. 

OBJECTIVES 

1)  Provide  satellite  image  analysis  and  meteorological  support  for  ONR  CDOM  cruises. 

2)  Provide  quantitative  data  synthesis  support  for  chemical,  biological,  and  physical  interactions  of 
CDOM  data. 

3)  Develop  ecological  equations  for  CDOM  cycling  in  the  coastal  marine  environment. 

APPROACH 

We  hypothesize  that  CDOM  cycling  is  a  deterministic  process,  one  that  can  be  explained  by  physical, 
chemical,  and  biological  interactions.  Furthermore,  coupling  experimental  data  with  environmental 
modeling  will  lead  to  the  development  of  a  set  of  ecological  equations  that  will  resolve  the  sources  and 
sinks  of  CDOM,  and  the  impacts  on  water  column  IOPs  and  AOPs.  A  previous  one-dimensional 
numerical  simulation  of  the  bio-optical  properties  of  the  Sargasso  Sea  (Ecological  Simulation  1 .0 
(EcoSim  1.0)  Bissett  et  ah,  1999a;  Bissett  et  ah,  1999b)  suggested  that  the  cycling  of  CDOM  could  be 
mathematically  described  and  validated.  In  this  case,  the  EcoSim  1 .0  results  were  validated  against  the 
multi-year  bio-optical  time  series  program  operating  at  the  Bermuda  Atlantic  Time-series  Station 
(BATS  Siegel  and  Michaels,  1996;  Siegel  et  ah,  1995).  In  particular,  simulated  CDOM  did  not  co-vary 
with  the  particulate  organic  absorption  or  the  chlorophyll  a  concentration,  as  was  previously  assumed 
in  these  oligotrophic  environments  (Gordon  and  Morel,  1983;  Smith  and  Baker,  1981). 

Separation  of  the  CDOM  signal  from  the  particulate  signal  in  oligotrophic  regions  is  important  for 
understanding  the  propagation  of  light  to  depth,  and  its  impact  on  water-leaving  radiance  signals  used 
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in  satellite  oceanography.  These  regions  cover  approximately  80%  of  the  world’s  oceans.  However,  the 
CDOM  signal  can  be  1-3  orders  of  magnitude  higher  on  coastal  shelves  and  near-shore  regions.  These 
higher  signals  result  in  part  from  the  higher  primary  productivity  found  in  coastal  regions,  but  also 
result  from  the  loading  of  CDOM  from  the  outflows  of  rivers  and  estuaries.  Prediction  of  in-water 
Inherent  Optical  Properties  [IOPs]  and  the  resultant  water-leaving  radiance  [Lw]  would  require  that  we 
accurately  address  the  ecological  cycling  of  CDOM,  i.e.  the  chemical,  photochemical,  biological,  and 
physical  sources  and  sinks  of  CDOM  in  the  near-shore  region. 

WORK  COMPLETED 

Mississippi  Cruise 

To  support  planning  for  the  CDOM  cruise  on  the  R/V  Pelican  of  April  5th- 1 1th,  we  began  data 
collection  on  March  28th.  Sources  for  the  data  included  the  processed  SeaWiFS  and  AVHRR  imagery 
from  NRL-Stennis  Code  7240  (http://www7240.nrlssc.navy.mil/),  COAMPS  27  km  current  and 
forecast  FNMOC  (http://152.80.49.205/PUBLIC/),  NCAR  Research  Application  Program 
(http://www.rap.ucar.edu/weather/),  NOAA  National  Weather  Service  (http://www.srh.noaa.gov/lix/), 
and  the  Colorado  Center  for  Astrodynamics  Research  Gulf  of  Mexico  Near  Real-Time  Altimeter  Data 
Geostrophic  Velocity  Viewer  (http://www-ccar.colorado.edu/~realtime/gom-real-time_vel/).  In 
reviewing  the  SeaWiFS  and  AVHRR  imagery  on  the  28th,  we  noticed  some  interesting  SST  and  color 
anomolies  with  definitive  edges  around  the  Mississippi  Delta  region  that  suggested  some  far- Held 
physical  forcing  may  be  impacting  the  plume  area  during  the  course  of  the  experiment.  Thus,  we 
proceeded  to  analyze  data  from  the  beginning  of  March  in  order  to  forecast  possible  directions  of  the 
Mississippi  plume  over  the  cruise  period. 

A  synthesis  of  the  satellite  and  weather  forecast  data  was  collected  each  day  into  PowerPoint  files 
following  the  28th  until  the  end  of  the  cruise  on  April  1 1th.  Each  day  these  syntheses  were  placed  on 
the  UMB  FTP  server  for  access  by  the  R/V  Pelican  scientific  personnel  and  other  interested  parties. 
These  syntheses  are  also  currently  available  on  request  from  FERI. 

EcoSim  2.0 

During  this  period  we  completed  a  re-analysis  of  the  color  cycling  equations  for  wavelengths  less  than 
440  mn  EcoSim  1.0.  These  equations  were  built  from  disparate  data  sets  of  DOM  and  CDOM  creation 
and  destruction  and  used  to  simulate  the  cycling  of  color  in  the  Sargasso  Sea.  A  critical  error  in  this 
simulation  was  the  extreme  loss  of  absorption  at  shorter  wavelengths  (<440  mu),  This  loss  was  thought 
to  result  from  a  number  of  possible  errors  in  the  quantitative  description  of  color  cycling  relating  to  -1) 
indirect  creation  of  color  by  phytoplankton  and  bacteria  lysis  and  grazing;  2)  direct  creation  of  color 
bacteria  excretion,  3)  destruction  of  color  via  heterotrophic  utilization  of  CDOM  for  energy 
requirements,  4)  destruction  of  color  via  photolysis  of  CDOM  to  colorless  DOM,  5)  destruction  of 
color  via  photolysis  of  CDOM  to  DIC. 

Each  of  the  above  factors  was  critically  analyzed  during  the  period  with  a  new  set  of  equations  and 
parameters  derived  for  the  CDOM  and  color  detritus  cycling.  These  new  equations  were  included  into 
a  larger  simulation  analysis  of  the  West  Florida  Shelf  (WFS)  during  the  fall  of  1998.  In  this  simulation 
work,  we  were  attempting  to  hindcast  color  and  red  tide  dynamics  observed  during  an 
EPA/NOAA/ONR/State  of  Florida  funded  field  campaign  (ECOHAB).  Multiple  simulations  were  run 
in  order  to  test  the  CDOM  equations  and  their  sensitivity  to  parameter  settings. 
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RESULTS 


Mississippi  Cruise 

Figure  1  shows  the  CCAR  geostrophic  flow  calculations  from  March  28th.  A  warm  core  ring/loop 
current  excursion  is  clearly  evident  to  the  south  of  the  delta.  The  strong  eastward  current  to  the  south 
of  the  delta  held  for  the  entire  cruise  period.  However,  the  plume  appears  to  be  more  responsive  to 
surface  winds  than  to  offshore  flows.  The  surface  winds  from  the  east-southeast  were  driven  by  a  low- 
pressure  center  to  the  north,  and  a  high-pressure  center  just  to  the  east,  which  for  the  most  part  kept  the 
plume  flowing  to  the  west-southwest  (Figure  2  and  3). 

We  also  found  that  the  SeaWiFS  backscattering  algorithm  at  555  nm  developed  by  B.  Arnone  of  NRL- 
Stennis  was  a  useful  indicator  of  surface  plume  dynamics  near  the  mouth  of  the  Mississippi  (Figure  3). 
In  some  cases  this  is  a  more  useful  algorithm  than  those  focused  on  other  SeaWiFS  bands.  This  may 
appear  to  be  somewhat  of  a  paradox,  given  that  the  CDOM  does  not  absorb  as  much  light  at  555  nm  as 
it  does  at  412  nm.  However,  the  Mississippi  River  also  has  a  high  suspended  solid  load,  which 
produces  a  high  Lw  signal  at  555  nm,  making  it  possible  to  readily  distinguish  the  plume  waters  from 
background  Gulf  of  Mexico  waters.  The  higher  absorption  of  light  at  412  nm  may  cause  signal-to- 
noise  difficulties  for  CDOM  retrieval  algorithms.  In  addition,  there  may  be  calibration  problems  with 
the  412  nm  channel  (per.  comm.  B.  Arnone).  For  this  river  plume  under  these  conditions,  it  may  be 
that  plume  dynamics  can  best  be  discerned  remotely  via  backscattering  calculations  versus  absorption 
calculations. 
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Figure  1.  TOPEXZERS2  analysis  of  sea  surface  height  and  calculated 
geostrophic  flow  on  March  28,  2001. 
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Figure  2.  FNMOC  COAMPS  27  km  resolution  surface  wind  forecast  for  April  5,  2001.  Winds 
around  the  Mississippi  Delta  are  from  the  east-southeast. 


Figure  3.  NRL  SeaWiFS  data  products  for  the  Mississippi  Bight  on  April  5,  2001.  Upper  figure  is 
adg(412 )  calculated  using  K.  Carder  (USF)  algorithm.  Lower  figure  is  the  bb(555)  calculated  using 

B.  Arnone  (NRL)  algorithm. 
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EcoSim  2.0 


Re-analysis  of  the  indirect  production  of  color  by  phytoplankton  and  bacteria  led  to  a  significantly 
reduced  value  of  the  fraction  of  CDOM  released  from  lysis  and  grazing.  The  fraction  of  carbon  was 
released  as  CDOM  was  reduced  from  5.5  to  3.23%  during  these  processes.  The  value  was  calculated 
during  a  reconstruction  of  the  labile  and  recalitrant  specific  absorption  (Bissett  et  al.,  1999a)  using  the 
specific  absorption  of  pure  water  (Pope  and  Fry,  1997)  and  the  assumption  of  net  non-use  of  color  for 
bacterial  heterotrophic  growth  (Moran  et  al.,  2000;  Vodacek  et  al.,  1997).  This  change  led  to  a 
significant  increase  in  the  quantity  (-300%  increase)  and  total  absorption  of  light  by  labile  CDOM. 
Prior  to  this  change  heterotrophic  bacteria  had  consumed  almost  all  of  the  labile  CDOM  production, 
because  in  oligotrophic  waters  they  are  typically  energy  (labile  DOM)  limited  (Bissett  et  al.,  1999b). 

The  other  significant  change  was  a  decrease  in  the  rates  of  photolysis.  The  initial  rates  (Bissett  et  al., 
1999a)  were  found  to  have  been  derived  with  experimental  data  that  was  not  appropriate  for  the  type  of 
calculations  we  were  attempting.  In  one  study,  it  appeared  that  there  were  biological  contaminants  in 
the  photolysis  experiment.  In  another,  the  decay  rate  was  calculated  over  the  entire  time  integral  (>8 
hours).  While  this  was  accurate  for  that  time  integral,  the  EcoSim  calculations  use  a  much  smaller  delta 
t  (<300  seconds),  and  as  a  result  our  estimated  decay  rate  was  much  too  fast  for  our  method  of 
calculation,  i.e.,  we  used  a  linearly  calculated  time  decay  in  an  exponential  type  equation.  Our 
photolysis  decay  rates  for  DIG  and  colorless  DOC  were  thus  reduced  (at  410  nm)  from  0.62  and  0.59 
|lM  m  hr'1,  respectively,  to  0.0193  and  0.0034  |iM  m  hr’1.  This  change  also  resulted  in  a  significant 
increase  in  the  accumulation  of  CDOM  (-70%)  and  an  increase  in  the  total  absorb tion  by  CDOM. 

These  new  terms  were  placed  into  our  2-D  simulations  of  the  West  Florida  Shelf  for  1998  (see 
progress  report  for  N00014-98-1-0844).  In  these  numerical  experiments  we  were  simulating  the 
impacts  of  esturine  discharges  of  33  psu  waters  onto  the  shelf  waters  of  36  psu.  Figure  4  shows  the 
total  absorption  for  these  at  412  nm.  These  compare  quite  well  with  the  AC-9  derived  total  absorption 
at  412  nm  versus  salinity  gradients  (Figure  5)  collected  during  the  ECOHAB  Process  cruises  (per. 
comm.  G.  Kirkpatrick,  Mote  Marine  Laboratory  and  O.  Schofield,  Rutgers  University) 

IMPACT/APPLICATIONS 

To  forecast  the  clarity  of  the  water  column  over  both  short-  and  long-tenn  time  horizons  requires  an 
accurate  quantification  of  the  ecological  cycling  of  CDOM.  Incorporation  of  a  validated  set  of  CDOM 
equations  into  a  larger  three-dimensional  ecological  simulation  will  increase  the  veracity  of  the 
predictions  of  Inherent  and  Apparent  Optical  Properties  [IOPs  and  AOPs],  and  help  achieve  the  goal  of 
forecasting  optical  properties  as  a  function  of  the  biological,  chemical,  and  physical  forcing. 

RELATED  PROJECTS 

1)  Robert  Chen  (UMB)  was  chief  scientist  on  the  ONR- funded  CDOM  cruise,  and  collects  the  spatial 
distribution  of  CDOM,  as  well  as  a  suite  of  other  environmental  variables  with  the  ECOShuttle. 

2)  Mary  Ann  Moran  (U  of  Georgia)  and  Richard  Zepp  (EPA)  are  conducting  laboratory  and  field 
measurements  on  the  photolysis  of  CDOM,  as  well  as  the  bacterial  utilization  of  in  situ  and 
photodegraded  CDOM. 
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3)  William  Miller  (Dalkousie,  CA)  is  conducting  experiments  on  the  photolysis  of  CDOM  as  a 
function  of  the  spectral  distribution  of  irradiance. 
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Figure  4.  Simulated  a/412)  from  2-D  West  Florida  Shelf  EcoSim  2.0  simulation  during  fall  1998. 
Cross  shelf  contours  show  impacts  of  low  salinity  waters  (33  psu)  on  total  absorption  resulting  from 

CDOM  influxes  at  shoreward  boundary. 


Figure  5.  ECOHAB  Process  Cruise  Data  from  1998  showing  salinity  versus 
absorption  relationship  from  AC-9  measurements  at  412  nm. 
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